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Abstract 

Stationary to pulsating soliton bifurcation analysis of the complex Ginzburg-Landau equation 
(CGLE) is presented. The analysis is based on a reduction from an infinite-dimensional dynamical 
dissipative system to a finite-dimensional model. Stationary solitons, with constant amplitude and 
width, are associated with fixed points in the model. For the first time, pulsating solitons are shown 
to be stable limit cycles in the finite-dimensional dynamical system. The boundaries between the 
two types of solutions are obtained approximately from the reduced model. These boundaries are 
reasonably close to those predicted by direct numerical simulations of the CGLE. 
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The complex Ginzburg-Landau equation (CGLE) is one of the basic equations for mod- 
elhng modulated amplitude waves [ij, spatio-temporal dynamics and spontaneous develop- 



ment of coherent structures in a variety of nonlinear dissipative systems 



Examples 



include pulse generation by passively mode-locked soliton lasers [4j, signal transmission in 
a..-opt,ca. commu„,catio„ Me. g travelling waves b,„a.,v flu.d ,„ixt„res Q. a„d also pat- 
tern formation in many other physical systems |3|- Complicated patterns consist of simpler 
localized solutions like fronts, pulses, sources and sinks jsf. 

Pulsating soliton solutions of dissipative systems have attracted a great deal of attention 
in recent years. They have been found numerically ^ IQ, ^ and observed experimentally 



111 ] in a fiber laser. Pulsating solitons form one set of possible localized solutions of the 
CGLE, and they exist on an equal basis with stationary solitons. Such localized waves exist, 
in various forms in biology, chemistry and physics. 

A pulsating soliton can be be described as a limit cycle of an infinite-dimensional dissi- 
pative dynamical system Q]. It is different from the higher-order solitons that are usually 
connected with an integrable model Although numerical simulations show clearly the 
existence of pulsating solutions and their bifurcations from stationary solitons, so far there 
has been no progress in finding analytic expressions for pulsating solutions and bifurca- 
tion boundaries. The problem is not simple as there are several parameters of the CGLE 
that define the regions of existence for both stationary and pulsating solitons. Hence, the 
bifurcation boundaries are surfaces in this multi-dimensional space of the parameters. 

In this work we use a reduction from an infinite-dimensional to a five-dimensional model, 
and we aim to find localized solutions of the CGLE and the transformations that they 
are subjected to when the system parameters are varied. Although exact solutions of the 
CGLE do exist jsl, they can be presented explicitly only for certain relations between the 
parameters of the equation. Furthermore, only stationary solutions can be found. Hence, 
we are faced with the necessity of finding an efficient approximation to tackle the problem. 
We have found that the method of moments, originally developed by Maimistov 14] for the 
perturbed nonlinear Schrodinger equation (NLSE) can be used for solving our problem. The 
moments are the integral characteristics of the field under consideration. In principle, there 
are an infinite number of equations for moments. One can obtain exact results by using the 
complete set of these equations. However, in practice, one uses a trial function with a finite 
number of parameters, and this is the way to obtain a significant reduction in the number 
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of variables used for the description of the dynamics. 

The cubic-quintic complex Ginzburg-Landau equation, in dimensionless form, is written 

as 

i6ij + ie\i)\^i) + i(3i)^^ + i/ilV-l V = RW\ > (1) 

where il){x,t) is the normalized envelope of the field, t and x are time and spatial variables, 
respectively, D is the group velocity dispersion coefficient, v is the parameter of the quintic 
nonlinearity, 5 represents the linear loss, e is the nonlinear gain coefficient, (3 stands for the 
spectral filtering, and /i characterizes the saturation of the nonlinear gain. Stable soliton 
solutions of the CGLE exist only for the following choices for the signs of the coefficients: 
5, /i < 0, /5, e > 0, and any sign for v and D (see e.g. 0, [l^). Hence, in this work, we limit 
ourselves only to this range. 

The method of moments is a reduction of the complete problem of the evolution of 
a field that has an infinite number of degrees of freedom to the evolution of a finite set of 
pulse characteristics. For a localized solution with a single maximum, these include the peak 
amplitude, pulse width and center-of-mass position. For an arbitrary localized field, the two 
integrals, namely the energy Q and momentum P 

/oo 1 /"^ 

i^pdx, p = - {^rx-r^x)dx. (2) 
-oo Z J— oo 

are two basic variables evolving with t. Three hi gher -order generalized moments, related to 



the pulse, are given by the following expressions 



h = x\ip\ dx , h = j (x - xo) dx , 

J —oo J —oo 

POO 

h= {x - Xo){'ilj*ip^ - ipipDdx , (3) 



where xo{t) = h/Q. Using Eq. one can derive the evolution equations for the generalized 
moments. For the five integrals given above, one can obtain the following jl^ : 

dQ 



—J- = i [ {'ipR* — ip*R)dx , 
dt J —oo 

/oo 
{i>xR* + ip*R)dx , 
-oo 



dP roo 

'dt " 



oo 

oo 



^ = iDP + il x{^R* -^*R)dx , (4) 
dt J —oo 



= -iDh + i {x- xofiipR* - ^*R)dx 

Cit J—OO 

dh 2P^ + z r {2Dm' - \^\')dx + 

Cit J—OO 



dt 
2i 



/oo roo 
{x - xo) {ipxR* + ^IR) + i / {i^R* + ^*R)dx . 
-oo J — OD 

Equations Q are general in the sense that they are vahd for a large class of NLSE-type 
evolution equations, including Eq. with arbitrary coefficients as a particular case. If we 
use an exact solution of Eq. (HJ, V') then equations (jH) are exact. 



De sufficient when one 
. In Ref. [isl, such an 



For problems of a certain class, even the first two Eqs. (@]) may 
deals with exact two-parameter reductions of the CGLE solutions [1 
approach was used to find the CGLE solutions in the form of stable soliton pairs and trains. 
The method of moments has also been applied to find stationary solutions of the CGLE 
m Ref. Ql 

although there was no attempt at finding pulsating solitons. Our tests showed 
that five is the minimum number of moments needed to describe pulsating solitons using 
the reduced system. Having more of them may improve the accuracy, but the complexity of 
the analysis then increases dramatically. 

The choice of the trial function is crucial for obtaining solutions with the desired prop- 
erties. Any reduction from an infinite-dimensional to a finite-dimensional system will have 
deficiencies. In approaches like this, the choice of the trial function can only be justified 
at the last stage of analysis, when the approximate solutions are compared with numerical 
simulations of the original equation. We use the fact that soliton solutions remain localized 
even when they are pulsating. Therefore, we take the sech-function: 

^{X, t)=A Sech f^—lR) gi[6(x-xo)+c{x-xo)2] ^ 



w 

where A{t), w{t) and xo{t) are the amplitude, width and maximum position of the pulse, 
respectively, b{t) is the soliton velocity and c(t) is the chirp parameter. The chirp is highly 
important, as the numerical simulations show. The number of parameters in the 

trial function must correspond to the number of moments used in the set of equations 
More complicated trial functions need more equations in 

Now, the generalized moments can be expressed in terms of the variable parameters of 
the trial function. Evaluation of the integrals Q and (jHl), with a help of Eq. gives the 
following expressions: 

Q = 2A^w, P = -2iA^wh, h = 2A'^wxo, 



(6) 



Then, using Eqs. (jl)), one can obtain a set of ordinary differential equations for the sohton 
parameters Q, w, c, Xq and b: 



Qt = F, = 2Q 



e Q 2/ig2 _ 5/5 



3 1(7 15w^ 
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Fixed points (FPs) of the dynamical system ((Zj) can be found from the set of algebraic 
equations = 0, j = 1, . . . , 5. The stability of the FPs is determined from the analysis of 
eigenvalues Aj, j = 1, . . . , 5, of the Jacobian matrix Mij = dFi/dpj, where {pi, . . . ,^5} = 
{Q, w, c, Xo, b}, and z = 1, . . . , 5. If the real part of at least one eigenvalue is positive, then 
the corresponding fixed point is unstable. In principle, the whole five-dimensional dynamical 
system ((Tj) can be studied using the specialized software 1^1 • We analyse the system based 
on general theory and using further simplifications. 

Firstly, we consider solutions with 6 = in Eqs. (|7j). The real part of the eigenvalue 
which corresponds to b is negative for any set of the system parameters. The soliton center, 
xo{t), takes a constant value for t —>■ 00 and the real part of the corresponding eigenvalue is 
zero (neutral stability). Therefore, we can make a further reduction and consider a system 
with only three variables, viz. Q,w, and c. We denote the three corresponding eigenvalues 
by Ai, A2 and A3. Note that, since the characteristic equation for these eigenvalues is cubic, 
then either Ai = A2 and A3 is real, or all three Aj are real. 

We find FPs numerically by solving the algebraic equations, and calculate the charac- 
teristic eigenvalues by analyzing the Jacobian matrix. In this way, we have identified the 
regions where solutions have distinctive features. The bifurcation diagram for the dynamical 



5 



system (|7j) in the {u, e)-plane is presented in Fig.[T] We take D = 1 and S = —0.1 henceforth 
in the paper. To facihtate comparison with exact results, the parameters of the model are 
chosen to be the same as in the numerical simulations of Ref. Q|, namely /x = —0.1 and 
p = 0.08. 

When the value of the gain e is small, there are no stable or unstable FPs (i. e. no 
stationary solitons) in the system. This region is located below e = 0.6, and therefore it is 
not shown in Fig. ^ For moderate values of e, there are two FPs. In the region below the 
curve 1 in Fig^, one FP is stable while the other one is unstable. The signs of the real 
and imaginary parts of the eigenvalues, for the stable FP in this region are the following: 
R-e[Ai 2] < 0, Im[Ai] = — Im[A2], Re[A3] < 0, and Im[A3] = 0. The stable FP corresponds to 
a stationary CGLE soliton with constant soliton parameters pj. 

Curve 1 in Fig^ is the bifurcation boundary (threshold) where the stable fixed point 
turns into an unstable one. Above curve 1, Re[Ai^2] > 0, Im[Ai] = — Im[A2], Re[A3] < and 
Im[A3] = 0. On line 1, a stable FP is transformed into a stable limit cycle (LC) and an 
unstable FP. We confirmed, numerically, that the stable limit cycle of model does indeed 
exist between the solid curves 1 and 2. Two examples of LC in (Q, w, c)-space are shown 
in Fig. I2K and ^p. The limit cycle has zero size on line 1 in Fig^. This corresponds to 
a super-critcal Hopf bifurcation The limit cycle elongates in Q axis direction when e 
increases. It becomes infinitely large and breaks off on line 2 in Fig^. 

The smaller region in Fig^ surrounded by the light solid curve, is copied from FigO^ 
for comparison. This is the boundary for the existence of pulsating solitons found by direct 
numerical simulations of CGLE (^. One can see that the thresholds (solid curves 1 and 2) 
of the existence of the stable LC provide reasonably good estimates for the boundaries. 

The period of oscillations, T, depends on the parameters of the system. A plot of T 
versus e at constant z/ = —0.09 is given in Figl^l The upper curve represents the results 
of exact numerical simulations of the CGLE, while the lower curve shows the period of 
oscillations found from our finite-dimensional model. There is an apparent difference in the 
values found for the period, due to the drastic reduction in the number of degrees of freedom 
in the model. However, the two curves have the same qualitative behaviour. In particular, 
each curve starts with a finite value of the period T at the lower boundary of the region 
where pulsating solitons exist, and increases to infinity at the upper boundary. 

Above the curve 2, the soliton amplitude A = [Q/(2w)Y^'^ remains almost constant. 
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but the soliton energy Q and the width w increase monotonically with t. This motion is 
tantamount to the locahzed solution of CGLE with constant amphtude and the width that 
increases indefinitely. The final stage of this motion is an asymptotic transformation of the 
soliton into two separating fronts with constant velocities. Thus, the qualitative agreement 
above the line 2 is also fairly good and the main features of the dynamics are captured 
correctly. 

Another slice of the parameter space, namely the {fi, e)-plane, is shown in Fig. As in 
Fig. there are two FPs in the region presented. One of them is always unstable. The 
stable FP exists below the solid curve 1. The bifurcation on the curve 1 is the same as 
that in Fig. The solid curve 2 corresponds to the transition between limit cycles and 
the solutions with increasing width. From Fig. ^ and Fig. EJ one can see that the results 
obtained from the model (|7j) for pulsating solitons and the boundaries for their existence are 
in reasonable agreement with direct numerical simulations of the CGLE over a wide range 
of the system parameters. 

Thus, considering the evolution of the dissipative soliton profile, the model ((Zj) is able to 
predict: (a) stable and unstable stationary solitons; (b) periodic soliton pulsations; (c) the 
unlimited increase in the soliton width at constant amplitude (which is equivalent to splitting 
solitons into moving fronts); and (d) bifurcations between these dynamical behaviours. 

In summary, we have derived a finite-dimensional dynamical system associated with the 
CGLE. We have demonstrated that the pulsating solitons of the CGLE correspond to the 
limit cycles of this system, while stationary solitons of the CGLE are related to the fixed 
points. We have found the approximate bifurcation boundaries between different types of 
solutions in the parameter space. 

This work was funded by the Australian Research Council. The authors are grateful to 
Dr. A. Ankiewicz for a critical reading of the manuscript. 
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FIG. 1: (a) Regions of existence and stability of FPs and limit cycles (LCs) of the reduced system 
in (z^, e)-plane. The central region, between the two solid lines 1 and 2, corresponds to stable LC. 
(b) Regions of existence for various soliton solutions obtained from numerical simulations of CGLE 
The region for pulsating solitons, found numerically in Ref. is copied from (b) to (a) for 
comparison. The parameters of the dynamical system are shown in (b). 




FIG. 2: Two examples of stable limit cycle of the system in (Q, c)-space. The parameters 
are shown in (a) and (b) respectively. Note the different scales for the axes in (a) and (b). 
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FIG. 3: Period of pulsations T versus e. Upper curve correspond to exact numerical results while 
the lower curve describes the data obtained from the low-dimensional model. 
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FIG. 4: (a) Diagram of existence and stability of FPs and LCs of the reduced model in the e)- 
plane. Stable LCs exist between the curves 1 and 2. (b) Existence diagram based on numerical 

n 

simulations of CGLE [10|]. The solid curve bounding the region for pulsating solitons in (b) is 
copied to (a) for comparison. The parameters of the dynamical system are shown in (b). 
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